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. Abstract 

| We investigate the fermionic one-loop cross section for the two photon annihilation of Kaluza- 

(N ■ 

Klein (KK) dark matter particles in a model of universal extra dimensions (UED). This process 

Ph. 

. gives a nearly mono-energetic gamma-ray line with energy equal to the KK dark matter particle 

I/") \ mass. We find that the cross section is large enough that if a continuum signature is detected, the 

energy distribution of gamma-rays should end at the particle mass with a peak that is visible for an 
energy resolution of the detector at the percent level. This would give an unmistakable signature 



of a dark matter origin of the gamma-rays, and a unique determination of the dark matter particle 
mass, which in the case studied should be around 800 GeV. Unlike the situation for supersymmetric 



■ models where the two-gamma peak may or may not be visible depending on parameters, this feature 



seems to be quite robust in UED models, and should be similar in other models where annihilation 



Q_i' into fermions is not helicity suppressed. The observability of the signal still depends on largely 

unknown astrophysical parameters related to the structure of the dark matter halo. If the dark 
matter near the galactic center is adiabatically contracted by the central star cluster, or if the dark 

X. 

matter halo has substructure surviving tidal effects, prospects for detection look promising. 
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I. INTRODUCTION 

The mystery of the dark matter, first indicated by Zwicky some 70 years ago Q], is still 
with us. Thanks to remarkably accurate measurements of the microwave background ^, 
supernova luminosity distances P] and the redshift distribution of baryonic structure [4] , we 
now have, however, a good handle on the required dark matter density of the Universe. 

Among the various possibilities for the dark matter constituents, weakly interacting mas- 
sive particles (WIMPs) are particularly attractive, as gauge couplings and masses in the 
electroweak symmetry breaking range (50 - 1000 GeV, say) give a thermal relic abundance 
of particles that naturally comes within the observed range. Supersymmetry or Kaluza- 
Klein (KK) particles in models with large extra dimensions are examples which span some 
of the range of the phenomenology to expect 

Dark matter particles can in principle be detected in three different ways: In accelerators, 
if the energy is high enough to produce them. Through direct detection in ultra-sensitive 
detectors, where particles from the local neighborhood of the halo deposit enough energy to 
be detected. Finally, indirectly through the contribution to the cosmic ray flux (positrons, 
antiprotons, neutrinos and photons of all energies up to the particle mass) caused by an- 
nihilations in the halo (or in celestial bodies for neutrinos). Here we will focus on indirect 
detection through high energy gamma-rays, a signal that has recently become much more 
interesting: because of new large-area air Cherenkov telescopes (ACTs) 0,3 and the launch 
of the Integral satellite ls|. In a couple of years, a new satellite, GLAST 3], will take the 
field of high energy gamma-ray astrophysics to yet another level of accuracy and precision. 

There are many reported indications of cosmic ray signals which could in principle at 
least partly be due to the annihilation products of dark matter particles. Examples are 
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possible excesses in microwave radiation 1()|. in nositrons [n, 12 _an< 
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gamma-rays from 



171 ] and TeV range 



the galactic center (GC) in the MeV range |l3l. Il4j. GeV range [15, 

[lsl . Indeed, there have even been suggestions that the combined data on positrons, 
gamma-rays and antiprotons could be better explained by adding a neutralino annihilation 
contribution j^J. 

The problem with these indications is that they lack distinctive features, which will make 
the alternative hypothesis that there are unknown background contributions difficult to 
reject. A spectacular signature would on the other hand be the mono-energetic gamma- 
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rays expected from WIMP annihilation into 77 |21| or Z7 |22j. The phenomen ology o f 
this process has been worked out in quite some detail for the supersymmetric case jl 71 . |23| . 
Here we present the first analytical and numerical calculations of the dominant fermion 
loop contributions to the process B^B^ — > 77 for KK dark matter. This dark matter 
particle has recently been shown to contain a useful experimental signature in the form of 



a gamma-ray spectrum 
drops at that energy 



;hat remains remarkably flat up to the KK mass, and then steeply 
24j . Here, we investigate whether there could also be a sharp peak at 
E 1 = M B m due to the two-gamma process. Needless to say, there is no known astrophysical 
background which would cause a similar structure, so the question is only if the particle 
physics annihilation cross section and the dark matter halo density are large enough to 
make the signal visible. 

The paper is organized as follows. We begin by introducing the model that gives rise 
to the KK dark matter candidate In Sections IIHI and IIVI we then discuss the basic 

properties of the cross section and give the analytical result for the fermionic contribution. 
We also comment on a full numerical evaluation of the process 

5(1)5(1) 

— > 77 and argue 

that the fermionic part gives the main contribution. In Section El we insert the analytical 
results into a model of the Milky Way halo that represents the current knowledge based on 
N-body simulations and baryonic contraction near the GC. We also make some speculations 
on the existence of dark matter substructure ( "clumps" ) and finally, Section IVIl contains 
our summary and conclusions. 

II. UNIVERSAL EXTRA DIMENSIONS 

The lightest KK particle (LKP) is the first viable particle dark matter candidate to arise 
from extra dimensional extensions of the standard model (SM) of particle physics. It appears 
in models of universal extra dimensions (UED) j^j, [3, 13 (f° r the ^ YS ^ proposal of TeV 
sized extra dimensions, see jssjl); where all SM fields propagate in the higher dimensional 
bulk, and is stable due to conserved KK parity, a remnant of KK mode number conservation. 
This is analogous to supersymmetric dark matter models, where conserved R-parity ensures 
the stability of, e.g., the neutralino. Contrary to the supersymmetric case however, the 
unknown parameter space is quite small and will be scanned throughout by next generation's 
accelerator experiments. 
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We will consider the simplest, five dimensional model with one UED compactified on an 
S l /Z2 orbifold of radius R. All SM fields are then accompanied by a tower of increasingly 
massive KK states and at tree-level, the nth KK mode mass is given by 

m(") = ^/(n/ J R)2 + m| w , (1) 

where toew is the corresponding zero mode mass. However, the identification of the LKP 
is nontrivial because radiative corrections to the mass spectrum of the first KK level are 
.ypically larger than the corresponding electroweak mass shifts. A one-loop calculation 
26| shows that the LKP is given by the first KK excitation of the photon, which is well 
approximated by the first KK mode of the hypercharge gauge boson B^ since the elec- 
troweak mixing angle for KK modes is effectively driven to zero. The relic density 
was determined in [221 . Depending on the exact form of the mass spectrum and the re- 
sulting coannihilation channels, the limit from the Wilkinson Microwave Anisotropy Probe 
(WMAP) Q of QcDwih 2 = 0.12 ± 0.02 corresponds to 0.5 TeV < m BW < 1 TeV. Here 
^cdm is the ratio of dark matter to critical density and h is the Hubble constant in units 
of 100 kms -1 Mpc -1 . Collider measurements of electroweak observables give a current con- 
strain, of > 0.3 TeV E Q, whereas LHC 8 hon,d p r obe jactitation ra dii up to 

1.5 Tev |ad). 

Although all SM fields come with a tower of KK states, only those which are even under 
Z 2 have a zero mode. Fermions in five dimensions are not chiral, but chiral zero modes can 
still be singled out by assigning different transformation properties to the SU{2) doublet 
and singlet, respectively. Denoting four dimensional spacetime coordinates by and the 
fifth coordinate by y, the fermion expansion in KK states is then given by 

M^,y) = ^=p^ (0) (^) + ^^k^ n H^)cos^ + p^(^)sin n,y 



where 



V^R T V 7 v^R^L v ' R ^ v > R 

v v n=l 

-=P^< V) + -7=5 E ^ B) M c °s ^ + PlWW) sin VI 

P _^7 5 



(3) 

Thus at a given KK level, the fermionic field content is doubled as compared to the SM. In 
the following we neglect the zero mode fermion mass, and set the radiative mass corrections 
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FIG. 1: Fermion box contributions to B^B^ — > 77 including the first level of KK excitations. 
Not shown are the additional nine diagrams that are obtained by crossing external momenta. 

to the singlet and doublet KK excitations equal. The mass eigenstates are then given by 

#) = ^W e } = -7 5 #, (4) 

which contributes a relative minus sign in the Feynman rules for S^^-fermion vertices. 

The number of scalar degrees of freedom is also doubled, since the fifth components of 
the gauge bosons transform as four dimensional scalars. These are required to be odd under 
Z2 by gauge invariance. The KK excitations of Z5 and Wf combine with those of the three 
Goldstone bosons of the SM to form three physical and three unphysical scalars, the latter 
giving masses to zj? and W^ n \ The KK photon acquires a mass by eating A$, while the 
Higgs excitations remain physical. 

III. THE flWBW 77 POLARIZATION TENSOR 

Let us now discuss the direct annihilation of particles into photons. Taking all 
momenta as ingoing, the amplitude for this process is given by 

M = ef (pOef (p 2 )ef (p 3 )ef (jpt) M mz ^{pi,p 2 ,Pz,Pi) , (5) 

where pi >2 and ^4 denote the four- momenta of the B^'s and photons with polarization 
vectors ei, .., 64, respectively. 

Let us first consider those contributions to the polarization tensor that do not contain 
scalars. At the one-loop level, and taking into account only the first level of KK-excitations, 
they are given by the fermion box diagrams that are shown in Fig. ^ Photons couple to the 
charge of the KK fermions by the usual vector coupling, while the couplings of the B^s to 
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fermions contains a vector as well as an axial vector part and is given by 

-^y^lVa + T 5 )^ +c.c. (6) 

and 

^fefVa-f^+cx., (7) 

respectively. Here, is a zero level fermion and ffl ) its first singlet (doublet) KK 
excitation with hypercharge Y s (Yd) as explained above. From charge conjugation invariance 
it follows that there is no total axial vector contribution to the process (i.e. no remaining 
7 5 in the trace). This means that the polarization tensor can be written as 

Ai" 1 " 2 " 3 " 4 = -ia em a Y Q 2 (Y* + Y%) (Mf 2 " 3 " 4 + M^ 2 ^ A ) , (8) 

where we have pulled out a common factor for later convenience, with a em = e 2 /4ir, ay = 
gy/^iT and Q being the charge of the fermions in the loop. In the above expression, A4 V 
describes the contributions from vector couplings only, whereas M. a comes from those terms 
that only contain B^ axial vector couplings. In our case, with the couplings given by © 
and ©, 

The typical velocity of WIMPs is about v ~ 10~ 3 , i.e. they are highly non-relativistic. 
The momenta of the incoming B^'s are therefore well approximated by p = p\ = p 2 = 
(m B {i),0). Taking into account momentum conservation, 

2pP+i%+i% = 0, (10) 

and transversality of the polarization vectors, 

€l ■ p = e 2 ■ p = e 3 ■ p 3 = e 4 ■ Pa = 0, (11) 

the Lorentz structure of the polarization tensor can be decomposed as 

m BW m BW m BW 

B B B 

m BW m BW m BW 

+—^- g^Ps'pT + Ci g^g^ 4 + C 2 f M j» w + C 3 g^g^s , (12) 

m B(i) 
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where A, B and C are dimensionless scalars that depend on the external momenta and the 
masses appearing in the loop. 

Bose symmetry prescribes that the tensor must be invariant under the permutations 
(pi,yUi) <-> (p 2 ,j« 2 ) and (P3, fa) <-> (P4, /%). This translates into 

B 2 = — B 4 , 
B3 = —B 5 , 

,4(3,4) = ,4(4,3), 

Bi(3,4) = Si(4,3), 

B 2 (3,4) = -53(4,3), 

5 6 (3,4) = 5 6 (4,3), 

C x (3,4) = d(4, 3) , 

C 2 (3,4) = C 3 (4,3), (13) 

where we have introduced the notation A(3, 4) = A(p,p 3 ,p4), ,4(4,3) = ,4(p,p4,p 3 ) etc. 

Because of gauge invariance the polarization tensor furthermore has to be transversal to 
the photon momenta, i.e. 

fc l fc 2 fc 4 Jy/l filfJ,2fi3Li4 — fc l fc 2 fc 3 F4 - /vl Ml/i2/i3M4 ~~ u - ^ lfJ 7 

This gives 
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A = _B 2 — -B4 — 2B, 
C2 — B 5 , 

C 3 = -S 3 . (15) 

With the help of Eqs. (fT3*|) and ()15j) one can express the whole tensor amplitude by 
means of only three form factors, which we choose to be Bi, B 2 and B§. The full analytical 
expressions for these functions are given in the Appendix . We actually calculated all the 
form factors A, B and C independently and used Eqs. (fT3*j) and ()15j) merely as a consistency 
check of our results. 
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FIG. 2: The annihilation rate for m B (i) = 0.8 TeV as a function of the mass shift between the 
and KK fermions. The dependence on the mass is simply given by av oc m 



IV. CROSS SECTION 



For low relative velocities t> between the incoming B^ 1 ' particles the annihilation cross 
section a diverges, but the annihilation rate is proportional to av and finite. Assuming the 
same mass shift for all KK fermions, we find after summing over final and averaging over 
initial states: 



av 



lU7rm 2 BW 



{3 |5x| 2 + 12 \B 2 \ 2 + 4 |S 6 | 2 - 4Re [B x (B* 2 + B* 6 )}} 



(16) 



where 



9eff 



Q 2 (n 2 + ^ 2 ) = y- (17) 
The sum in (fTTj) is over all standard model charged fermions. From ordinary running of 
couplings, one expects l/a em (l TeV) ~ 123 and l/ay(l TeV) ~ 95. In models with extra 
dimensions one generically expects power-law rather than logarithmic running [jyl, which 
would further enhance the cross-section. However, the deviation from the standard case 
only becomes significant at energy scales considerably higher than the first KK level, so we 
neglect this effect here. Fig.|2]shows the annihilation rate as a function of the mass difference 
between the B^ and the KK fermions. 

Until now, we have only considered diagrams containing fermion loops. There are consid- 
erably more diagrams if one also takes into account the scalar sector of the theory; neglecting 
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those that are related to each other by obvious symmetries, there are actually a total number 
of 22 different types of diagrams with first KK level modes. Using the FormCalc package 
(3^], we implemented all these contributions numerically and found that they merely make 
up a few percent of the total cross section. Comparing this with the direct annihilation of 



neutralinos XX ~^ 77 2Jj> this is actually the order of magnitude one would naively expect 
for the total cross section when taking into account that m x /m B (i) ~ 10~ 2 . The large con- 
tribution from the fermion loops thus seems to be a special feature of the KK model we have 
considered here. However, it is well known that neutralino annihilation into charged leptons 
is helicity suppressed and the process XX ~~ > 77 is n °t dominated by fermions, either. This 
indicates that the small fermion loop contribution is a peculiarity for supersymmetry rather 
than its absence a special feature of the KK model. 

Higher KK levels contribute with additional diagrams, but these are suppressed by the 
larger KK masses in the propagators. Numerically, we have checked that adding second 
KK level fermions only changes our previous results by a few percent. Another possibly 
important contribution comes from diagrams with second KK level Higgs bosons H^> in the 
s-channel. We found, however, that for realistic mass shifts one is too far away from the 
resonance for these processes to be important. In conclusion, the analytical expression 
(|16[) for the fermion box diagrams shown in Fig. ^ is a good approximation for the whole an- 
nihilation process 

5(1)5(1) — > 77 and we expect the effect of other diagrams to be negligible 
compared to the astrophysical uncertainties described in the next section. 

One should also note that there are two other final states of B^'B^ 1 ' annihilation, Z7 
and i?7, which serve as sources of monochromatic photons. Due to the high mass of the 
these can most likely not be discriminated from those from two-photon annihilation. 
The fermionic contribution to the B^B^ — > Z7 cross section is easy to compute for a 
given diagram. The only difference compared to the two photon case is that the Z coupling 
to fermions contains both a vector and an axial vector. From the values of these couplings 
we have estimated, and checked numerically, that this gives a 10% enhancement of the 
gamma-ray signal coming from 77 annihilation alone. The if 7 line is left for future studies. 
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V. OBSERVATIONAL PROSPECTS 



The detection rate of annihilations in the GC depends strongly on the details of 
the galactic halo profile which, unfortunately, is to a large extent unknown. High resolution 
Nl-body simulations favor cuspy halos, with radial density distributions ranging from r _1 
to r -L5 [34] in the inner regions. It has further been suggested that adiabatic accretion 
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of dark matter onto the massive black hole in the center of the Milky Way may produce a 



dense spike of r 2 4 |35j]. This has, however, been contested |36|,|37|. On the other hand, 
adiabatic contraction from the dense stellar cluster, which is measured to exist near the GC, 
is likely to _ the halo dark matter substantially HQ- This ,„ea n8 that there is 
support for a rather dense halo profile very near the center - something that may be tested 
with the new 30m-class telescopes Bearing these uncertainties in mind, we will use the 
moderately cuspy (V" 1 ) profile by Navarro, Frenk and White (NFW) 



33 



The gamma-ray flux from WIMP annihilation in the GC is given by 



$ 7 (A0) ~ 2.92 • 10 



-n 



av 



10" 29 cm 3 s" 



\ x ^^vy (jGc)AnAOm -2 S -i ; (18 

V V m BW J 



where av is the annihilation rate, and (Jgc) ad. is a dimensionless line-of-sight integral av- 
eraged over AO, the angular acceptance of the detector. This way of writing the flux is 
convenient because it separates the particle physics of the particular WIMP candidate from 
the details of t he g alactic halo profile. For a NFW profile, one expects (Jgc) An AO = 0.13 6 



for AO 
follow 
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4l|, which is appropriate for example for the H.E.S.S. telescope (7|. Here, we 



24j and allow for an explicit boost factor b that parametrizes deviations from a pure 
NFW profile (b = 1) and may be as high as 1000 when taking into account the expected 



effects of adiabatic compression 

The direct annihilation B^B^ — > 77 produces mono-energetic photons with an energy 
that is to a very good approximation given by Ej = E B (v> ~ Tn B (i). Due to the velocities of 
the B' 1 ' particles, however, the emitted photons become Doppler-shifted and the observed 
signal acquires a natural linewidth of v ~ 10~ 3 . Current ACTs have energy resolutions of 
about 10 to 20 percent . Satellite-bourne telescopes have better spectral resolution but 



are limited in energy range; neither Integral (E^ < 10 MeV) [8( nor GLAST (E 1 < 300 GeV) 
[t| can reach the energies of some hundred GeV to a few TeV that are required in order 
to study the expected signatures of KK dark matter. The next generation of high energy 
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FIG. 3: The continuous gamma-ray flux that is expected from KK dark matter |2J] is plotted as a 
solid line. To this, we have added the flux from direct annihilation as seen by a detector with an 
energy resolution of 2a = 2 % (dashed), 1 % (dotted) and 0.5 % (dash-dotted), respectively. The 
actual linewidth of the signal is about 10 -3 , with a peak value of 1.5 • 10~ 7 m™ 2 s~ x TeV -1 . For all 
cases, we have assumed m B (i) = 800 GeV, a mass shift m^i) /m B (i) = 1.05, and a moderate boost 
factor of b = 100. 

gamma-ray telescopes can be expected to reach even better energy resolutions, also in the 
TeV range, something which opens up exciting possibilities for the study of line signals from 



dark matter annihilations as the one described here 



21. 



Recently, the continuous gamma-ray spectrum from KK dark matter annihilation was 



found to have a very promising signature 



24j . an almost flat energy distribution that 



abruptly falls off at m B (i). Including the line signal that we have studied here, one finds a 
total spectrum that has an even more distinctive signature. This is shown in Fig. El As one 
can see, the next generation of experiments may well be able to see it - or falsify the KK 
dark matter hypothesis. Note that the relative strength of the line signal as compared to 
the continuous signal does not depend on the astrophysical uncertainties encoded in b. 

In the particular model we are looking at, a dark matter particle mass as low as 300 GeV 
is disfavored by the WMAP relic density bound. However, it could be that the does 
not make up all of the dark matter; one could also imagine a different model, with equally 
non-suppressed couplings to fermions, that predicts a somewhat higher relic density. In both 
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cases one might be able to push the dark matter particle mass down to a value such that 
an annihilation signal like in Fig. |3] could be seen already with the GLAST satellite, taking 
into account its projected spectral resolution at the highest energies. 



VI. CONCLUSIONS 

We have studied the possible distinct signature of annihilation into mono-energetic 
gamma-rays of KK dark matter, one of the most interesting dark matter candidates that 
has appeared in recent years. We have found that with a detector of percent level energy 
resolution, something that is difficult but not impossible to achieve, such a "smoking gun" 
signal could indeed be detectable. The mass range of the KK models that derive from UED, 
around 0.5 - 1 TeV, places the new generation of ACTs at the center of attention. If the 
halo density profile has a slope of 1/r or more and is furthermore processed by interaction 
with the baryonic matter near the center of the galaxy, then this region is the obvious place 
to search for the gamma-ray line. On the other hand, it may be that the galaxy halo is 
clumpy 42I and then the best choice may be to aim at the nearest dark matter clump. The 
upcoming satellite experiment GLAST will have difficulties detecting energies above 300 
GeV, but may do an all-sky search for the continuum signal expected in these models. The 
corresponding point sources without optical counterpart would be obvious targets for the 
observing programs of the different ACTs then in operation. An observation of the photon 
line would then unambiguously identify the dark matter origin, as well as give the mass of 
the lightest KK particle. 

Much of the phenomenology would be the same for other models which do not involve 
Majorana fermions as dark matter. That exception happens on the other hand for neutrali- 
nos in supersymmetric models. As the Majorana nature means that the annihilation into 
fermions is helicity suppressed, one has to rely on a strong coupling to W ± in the loop to 
obtain large rates for the gamma line. Although there are examples in supersymmetric mod- 
els of such a coupling, for example in the pure higgsino limit [12J or in split supersymmetry 
it is not a generic feature. 

In models with non-suppressed couplings to fermions, on the other hand, the calculation 
of the continuous gamma-ray spectrum as given in j^J would be essentially unchanged, as 
would the present loop calculation of the gamma line rate. In particular, the relative weight 
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of the two processes should be rather universal, meaning that whenever a continuous signal 
is claimed to have been detected, there has to be a bump at the dark matter particle mass 
which should stand out if the energy resolution is a few percent, or better. 

Finally, one should note that the annihilation channels B^'B 1 - 1 ' — > Z'j and B^B^ — > 
H'j has a similar structure to the one we have studied here and is another source of monochro- 
matic photons. We have estimated the Z7 process to give at least a 10% enhancement of 
the signal that we have discussed in this article. We leave the interesting issue of the Hj 
line open for later studies. 
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Appendix 

In the calculation of the loop diagrams, the following types of four-point functions appear: 
D ; D^; D^; D^^ (k 1} fc 2 , h; m 1} m 2 , m 3 , m 4 ) 

f d^g 1; g^, g^; g^g P ] <yhAtfh (iq) 

J nr 2 [g 2 - mj] [(q + hf - m|] [(q + k x + k 2 f - m|] \{q + k 1 + k 2 + k 3 ) 2 - m|] ' 1 } 

As is well-known, all these tensor integrals Dij can be reduced to scalar loop integrals [44], 
for which closed expressions exist Q|. In certain kinematical regimes, however, the original 
reduction scheme 44| of Passarino and Veltman breaks down. This is for example the 
case for the process that we are interested in, where the two incoming B^ particles have 
essentially identical momenta. We use the algebraic reduction scheme implemented in the 
LERG program that adopts an extended Passarino- Veltman scheme which can cope 
with special situations like this. 



We find that all the integrals one has to calculate for the first KK level 
can be reduced to only six basic scalar integrals. Using the conventions from 



ermion loops 



45 



46j (note, 



however, that our signature is different from theirs), the form factors defined in section ITTT1 
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become 

B l = -^-C (4, 0, 0; 77, 77, 77) - -^C (l, 0, -1; 0, 77, 77) - ff (4; 0, 0) 

77 - 1 7/ - 1 3(77+1) 

' t^j^ ") + 2 ' 5 '^^- 19 » B „ (1 ; 0. ,) - Mb(-l; 0, ,) + § , (20) 
B 2 = 4 '' ( '' + , 2) Co(4, 0, 0; n , n ) - 4 '' (2 '' + 1} C„(l, 0, -1; 0, r,, n ) - — -!-£>,, (4; 0, 0) 

1] — 1 7/ — 1 3(7/ + 1) 

+!fi? *<«™> - 13 'U^ 13 bm + BM) - 1 • (21) 
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MV + Z)^^ nn v , 4(377 2 + 77 -n , t n 1>n , 8 



^-^C (4, 0, 0; 77, 77, 77) + ^C (l, 0, -1; 0, 77, 77) - - B (A; 0, 0) 

77 — 1 77 — 1 3(77 + 1) 

2 < 13 » - 4 W„) + 41 "L +26 " 1 : 31 ^(1;0,„) - 5ft<-l;0,„) + I, (22) 



3(77-1) - U ^''"'" 1 3(77 2 -l) — 1 3 

where 77 = |^ m J — j . We calculated the scalar loop integrals appearing in the expressions 
above as follows: 

C (4,0,0;77, 77, 77) = ~ arctan 2 (^fi=f) > ( 23 ) 

C (l, 0, -1; 0, 77, 77) = i |Li 2 (-^j - Li 2 j , (24) 

5 (4;0,0) = 2-21og2 + i7r, (25) 

5 (4;77, 77) = 2 - log 77 - 2^7? - 1 arctan ( ^=f ) > ( 26 ) 

50(150,77) = 2-77log77 +(77 -l)log(77-l), (27) 

£?o(-l;0,77) = 2 + 77 log 77 -(77 + l)log(77 + l), (28) 

where the Spence function or dilogarithm is defined by 



i^-jf d.lS5e^ = f^. (29) 
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